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, Abstract. We study phase transitions of a system of particles on the one-dimensional integer 
lattice moving with constant acceleration, with a collision law respecting slower particles. This 
jyy • simple deterministic "particle-hopping" traffic flow model being a straightforward generalization 
Q . to the well known Nagel-Schreckenberg model covers also a more recent slow-to-start model as 
(-H ! a special case. The model has two distinct ergodic (unmixed) phases with two critical values. 

When traffic density is below the lowest critical value, the steady state of the model corresponds 
to the "free-flowing" (or "gaseous") phase. When the density exceeds the second critical value the 
model produces large, persistent, well-defined traffic jams, which correspond to the "jammed" 
^ (or "liquid") phase. Between the two critical values each of these phases may take place, which 

^ ' can be interpreted as an "overcooled gas" phase when a small perturbation can change drastically 
CD . gas into liquid. Mathematical analysis is accomplished in part by the exact derivation of the 
I life-time of individual traffic jams for a given configuration of particles. 

OO , 
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1 Introduction 



By a traffic flow we shall mean a collection of particles moving along a straight line according to 
^ . their velocities, and the law describing how those velocities are changing (i.e., the acceleration or 
^ I deceleration of particles) is called the traffic flow model. For a long time theoretical analysis of 
traffic flow phenomena has been dominated by hydrodynamic models in analogy to the dynamics 
of viscous fluid (see, e.g. [HI El EI and further references therein). The main problem with this 
approach is that it gives no information about the behavior of individual vehicles. Moreover, 
the discrete nature of the traffic flow has some features (like traffic jams) which do not have 
exact counterparts in the hydrodynamic interpretation. Additionally, considered as a kind of 
viscous fluid, the traffic flow turns out to have very peculiar properties. First, it is very much 
'compressible': the distance between two consecutive particles being originally quite far from 
each other may shrink to zero under the dynamics due to a large jam ahead of them. On 
the other hand, the opposite property does not hold: the distance between two consecutive 
particles can be enlarged at most by a constant depending on the parameters of the model but 
not on time. Mention also that those hydrodynamic models are difficult to treat in computer 
simulations of large networks, while it is hard to compare parameters of the models with empirical 
investigations. To overcome these difficulties cellular automata models have been invented in 
the beginning of 90s (see HJISj for reviews). 
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Figure 1: Freely accelerating dynamics of particles. i± and i'j^ denote the positions of neighboring 
particles at moments of time t and t + l and x± ~ corresponding velocities at time t. 



Consider a system of particles on the integer lattice Z moving with constant acceleration 
o € in a discrete time with a collision law respecting slower particles. To make this description 
precise we need a few definitions. Each particle is described by the pair where i G 

represents the position of the particle and Xi € its velocity, which might be both positive 
and negative. Fix a configuration of particles. For each i G by i_ < i we denote the site 
containing the particle from our configuration closest to the site i from the left side, and by 

> i the one containing the particle closest to i from the right side (see Fig. Note the 
asymmetry of the definition: i_ might be equal to i, while i+ is strictly larger than i. 

We shall say that a configuration of particles is admissible if for any i € Z"*^ we have: 

i+ - > max{0, Xj., - Xi^, Xi_ - Xi^], |a;i±|<v, (1.1) 

where v > is a parameter of the system describing the maximal allowed velocity. In other 
words all particles in an admissible configuration can be moved in arbitrary order by the distance 
equal to the corresponding velocities without collisions. 

On admissible configurations the dynamics is defined as follows. First we change simultane- 
ously the coordinates (i,Xj) of all particles according to the rules: 

i -^i+\xi\ (1.2) 
Xi ^min{xj + o, v} (1-3) 

where by [-J we denote the integer part of a number, i.e. [zj := max{n E Z : n < z}. Denote 
also by [•] the smallest integer not smaller than the considered number, i.e. [z] := min{n G 
Z : n>z}. 

After this operation for each particle (i, Xi) for which the admissibility condition (jl.lj) breaks 
down (which correspond to a 'collision' at the next step of the dynamics) we correct its velocity: 

Xi —I- minjxj, i+ — — 1}. (1-4) 

Then for all particles where the velocities remain negative we make an additional correction: 

Xi —f max{xj, i — {i — 1)_ — 1 — max{0, X(j„i)_ }}. (1-5) 

It is straightforward to check that after these corrections the configuration of particles becomes 
admissible.^ 

Fig. 121 illustrates the dynamics of our model in the case a = 1/2 and v = 1 . Here and in the 
sequel we mark the positions of particles in the configuration by their velocities and the positions 
of holes by dots. The first line in the figure corresponds to a finite segment of length 21 in the 
initial configuration (i.e. at time 0) and the subsequent lines describe the same segment for the 
moments of time from 1 to 6. Note that we have chosen the segment such that the particles not 

^We emphcisize that the update rules 11.211 - 11.51 are appUed simultaneously to all particles in a given 
configuration. 
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Figure 2: An example of the dynamics with a = ^, v = 1. The positions of particles are marked 
by their velocities and the positions of holes by dots. 



shown in the figure do not influence the dynamics of the ones in the segment under consideration 
during the first six time steps. 

Due to the constant acceleration after a finite time (of order v/a) all particles will start 
moving in the same (positive) direction. Therefore, since we shall be interested mainly in the 
limit behavior of the system as time goes to infinity, to simplify notation we shall assume from 
now on that velocities of all particles are nonnegative^ and shall describe the system by a 
configuration x = {xi }i G , where a nonnegative entry Xi corresponds to a particle at the site 
i with the velocity Xi > 0, while all other sites (containing 'holes') are set to —1 (having in mind 
that the 'holes' are moving in the opposite direction to the particles). Under this assumption a 
configuration x is admissible if for any i £ we have 

i-+Xi_<i+, Xi G [0,v] U {-1}. (1.6) 

The space of admissible configurations we denote by X. Using this notation the dynamics can 
be defined as a superposition of two maps T := A o a, where the maps a : X — > M°° and 
A : M°° — > X are defined as follows: 

■= 1-1 otherwise ' ^^'^^ 

(A^). _ fmin{xi + a, i+ - i - 1, v} if Xj > 
\ — l otherwise 

The map a is a (nonuniform) shift map, describing the simultaneous free shift of all particles 
in the configuration x by distances equal to their velocities, while the second map A describes 
the process of the acceleration/deceleration of particles (see Fig. |2|. Observe that after the 
application of the map A any configuration becomes admissible which ensures that the map T 
preserves the space of admissible configurations X. 

The dynamics described by the map T for integer values of the parameter a > coincides 
with the deterministic version of the well known Nagel-Schreckenberg model (see e.g., jlHH] for 
reviews). Note also that a version of the Nagel-Schreckenberg model for continuous values of 
the velocities and spatial coordinates was discussed in 0- 

We allow the velocity of a particle to be any real number between and v, but since the 
particles are moving on the integer lattice the actual shift of a particle (described by the map a) 
is an integer, which means that a particle starts moving only when its velocity becomes greater 
or equal to 1. 

One of the most striking observations in the theory of traffic flows is the fact that in a 
large number of models as time goes to infinity the limit 'average velocity' of the flow depends 
only on its density and does not depend on the finer characteristics. This dependence is called 

^Note, however, that due to the argument above all results in the paper remain valid in the initial setting with 
particle velocities of both signs. 
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Figure 3: Fundamental diagram for T: dependence of the limit average velocity V on the density 
of particles p. Here w := [l/a]- 

'fundamental diagram' in physical literature. Originally this observation has been established 
numerically for finite lattices with periodic boundary conditions, or in our terms for space pe- 
riodic configurations.^ Later this result has been proven rigorously for a broad class of models 
[HOIE]) assuming that the 'average velocity' is, indeed, the average velocity along a configu- 
ration, which we shall call the 'space average velocity'. 

Recently a number of nonmarkovian traffic flow models demonstrating the so called 'metastable 
states' (see, e.g. ^1 |Hj and references therein) have been introduced. Numerous numerical 
simulations of those models show that in distinction to the previous models the 'space average 
velocity' wildly fluctuates in time. Our own numerical results also show the same phenomenon 
with fluctuations up to 100% of the 'average' value if the acceleration a < 1. To demonstrate 
this choose a space periodic configuration with the spatial period of length 6 consisting of the 
periodically repeating pattern [01.. 1.]. Considering the dynamics with a = ^, v = 1 of only 
the main period of the configuration we get: 01. .1. — > a.l..O — > 1..1.0 — > .1..0a As we 

see starting from the 2nd iterate the 'space average velocity' fluctuates in time periodically. In 
general the limit as time goes to infinity 'space average velocity' may not exist at all, moreover, 
both upper and lower limits may fluctuate in time as well. This shows that such approach is not 
adequate here. On the other hand, the 'space average velocity' tells almost nothing about the 
movement of individual particles. Therefore it looks reasonable to consider the averaging in time 
for individual particles instead of the averaging in space. Note that due to the nonergodicity of 
the process under study these two quantities need not to coincide. On the other hand, it will be 
shown that the 'time average velocity' behaves much better and the main result of this paper 
is the proof that the limit 'time average velocity' V{x) for an admissible regular^ configuration 
X ^ X having a density p is described by the following multivalued function 

FD^v(p) = min{v, "J^^} U l[(i+ri/alv)-i, (i+v)-i] (z^) • v, (1.9) 

corresponding to the fundamental diagram shown in Fig. |31 The last term in the above relation 
describes the upper (unstable) branch of the fundamental diagram. Details related to lower 
and upper densities/velocities and non-regular configurations will be discussed in Sectional 
Let me emphasize that non-regular configurations are not excluded and results related to the 
fundamental diagram are obtained for all configurations having densities. 

Let us mention connections between the model under study and some previously considered 
cases. In the case when all velocities of particles are positive integers and the acceleration a is 
an integer as well we immediately recover the deterministic version of the well known Nagel- 
Schreckenberg model introduced in ^M,- If one assumes that a > v then our system coincides 

^Indeed, an admissible configuration with a spatial period i remains periodic in space with the same spatial 
period I which explains the connection to numerical results. 
''See definition in Sectional 



with the one studied analytically in [21 Ej . It is known that in the case a > v the asymptotic 
behavior is much simpler, in particular, only free flowing particles or free flowing holes may 
show up, and, moreover, traffic jams cannot grow in time. There is also another more technical 
difference due to the fact that in that case the dynamics of holes is completely symmetric to 
the dynamics of particles. This symmetry was heavily used in |21 01 for the analysis of the 
dynamics of high density particle configurations. If a < 1 there is no independent dynamics of 
holes as such: to define their dynamics one inevitably needs to take into account the velocities 
of particles. Therefore the analysis of a very complicated dynamics of high density particle 
configurations cannot be reduced to the low density case. Fortunately, as we shall see in the 
'jammed' phase the dynamics of individual holes becomes periodic in time (but not in space) 
after some initial stage. We shall see that if a < 1 then when the particle density exceeds 
some critical value the traffic jams can grow in time and, moreover, an arbitrary small jam (say, 
consisting of a few particles) may become arbitrary large with time. 

It is worth note that when the reciprocal to the acceleration a is an integer and the initial 
configuration has only integer velocities the model we are discussing coincides with the recent 
'slow-to-start' model, which was extensively studied on the numerical level (see e.g. jlfij, ^ 
and references therein), but no analytical results have been obtained so far. This model (being 
non-markovian in the sense that the dynamics depends on the history) turned out to be the first 
one demonstrating large fluctuations in 'space average' statistics or 'metastable states'. 

In this paper we restrict the analysis to the case of slow particles with v = 1, because the 
analysis of the fast particles with v > 1 and/or several lanes uses a rather different mathematical 
apparatus - substitution dynamics. Moreover, in the general fast particles case (when v > 
1) additionally to (static) traffic jams that we consider in this paper there are new types of 
(dynamic) jams when all particles in a jam might have positive velocities and are moving as a 
whole with the velocity strictly less than v. These new types of jams lead to additional branches 
in the fundamental diagram and will be discussed elsewhere. Still since some of technical results 
obtained in the paper remain valid for any v € we shall keep the notation v throughout the 
paper and shall specify explicitly if only the case v = 1 is considered. 

The paper is organized as follows. In the next Section we introduce basic notions including 
the notion of the traffic jam and its characteristics such as a basin of attraction and a life-time, 
and derive results describing their exact values as functions of a given configuration. Using these 
results in Section O we prove the validity of the fundamental diagram and provide its stability 
analysis in Section jll An alternative (more intuitive) proof of the validity of the fundamental 
diagram via space-time averaging is given in Appendix. 

2 Dynamics of traffic jams 

For a finite segment x[n,m], n < m £ Z,\_ U {0} =: Z>q of a configuration x G X we define the 
density of particles p{x[n,m]) as the number of particles in the segment Np{x[n,m]) divided by 
the length of the segment, i.e. 

pxn,m := ^ ' \' = Y,^>oix^), 2.1 

m — n + I m — n + I ^-^ ~ 

i=n 

where l>o(') is the indicator function of the set of nonnegative numbers. The generalization 
of this notion for the entire infinite configuration x £ X leads to the notion of lower/upper 
densities: 

p±{x):=lim(^'J^p{x[l,n]), (2.2) 

where (and in the sequel) limsup corresponds to the index '-|-' and liminf to the index '-'. 
If the lower and upper densities coincide their common value p(-) will be called the density 



of the configuration. According to Birkhoff ergodic theorem for any translationally invariant 
measure on the integer lattice the set of configurations having densities is the set of full measure. 
In particular, for any space periodic configuration the density is well defined.^ Probably the 
analysis of configurations with well defined densities would suffice, but since some configurations 
with differing lower and upper densities lead to interesting behavior we shall treat them as well. 

The reason to fix m = 1 in the definition above is that for any given integer m the limit 
points as n — > oo of the sequence {p{x[m,n])}n are the same and do not depend on m. The 
asymmetry with respect to the left and right 'tails' of the configuration reflects the fact that 
roles of those 'tails' are rather different in the long term dynamics. In particular, consider two 
symmetric (to each other) configurations: y^^^ having holes at all negative sites and particles 
at all positive ones, and y^^^ having particles at all negative sites and holes at all positive ones. 
Any particle in the configuration y^^^ has zero velocity and will never move; on the other hand, 
any given particle in y^^^ will eventually with time get the largest possible velocity v. 

Now we are ready to define time and space average velocities. Since the case when there are 
no particles in the configurations does not make any sense from the dynamical point of view, 
from now on we shall assume that any admissible configuration contain at least a finite number 
of particles. Let L{x, i, t) be the distance covered during t time steps by the particle initially 
(at t = 0) located at the site i_ of the configuration x. Denote 

V^'^^^x, i, t) = Vix, i, t) := iL(x, i, t) (2.3) 

V±{x,i) := hm ^^J" ^V{x,i,t) (2.4) 

^ m 

F^P-^(x[-n,m]) := -— ^ ^ ^ l>o(x,) • [x,\ (2.5) 

V^^^\x) := lim ( ^""P W^P^™(x[-n,m]). (2.6) 

n,in—*QO \ mf J 

In other words V^p^^^ describes the average velocity along all particles in a given configura- 
tion, while V corresponds to the average in time velocity of a given particle. More precisely 
these quantities describe average effective velocities since we take into account only the actual 
movement of particles and thus drop fractional parts of their velocities. 

As we shall show the time average statistics (in distinction to the space average one), does 
converge to the limit as time goes to infinity and the result does not depend on the initial site 
i, i.e. V{x, i) = V{x). 

To this end, let us show that the upper and lower densities are invariant with respect to 
dynamics. 



Lemma 2.1 p±{Tx) = p±{x) for any configuration x (z X. 

Proof. For any n,m (z Z>o we have 

\Npix[-n,m]) - Np{{Tx)[-n,m])\ < 2, (2.7) 

since during one iteration of the map at most one particle can enter the interval of sites from 
— n to m (from behind) and at most one particle can leave this interval.^ 

^Thus in numerical simulations only configurations with well defined densities can be observed. 
®If only nonnegative velocities are considered 1 (instead of 2) would sufiice in the inequality l|2.7^ , but in order 
to be valid in the original setting with velocities of both signs we put 2 here. 



3 — ^oo 

By the definition of the lower density there is a sequence of integers nj — > oo such that 

Np{x[l,nj]) j^oo 
rij 

On the other hand, due to the inequahty H2.7|) it follows that P-{x, ) is a limit point for partial 
sums for the configuration Tx. Therefore we need to show only that this is, indeed, the lower 
limit. Assume, on the contrary, that there is another limit point, call it for the partial sums 
for Tx such that ^ < /j_(x). Doing the same operations with the partial sums for Tx converging 
to ^ we can show that this value is also a limit point for the partial sums for the configuration 
x, and, hence, ^ cannot be smaller than p-{x). 

The proof for the upper density follows from the similar argument. □ 



Introduce a weight function: 

[U otherwise 

Referring to this function we shall drop the dependence on z if z = to simplify notation, i.e. 
w{d) = w. 

Lemma 2.2 Let x E X and let the limit V{x,i) be well defined. Then for any j G the limit 
V{x^j) is also well defined and coincides with V{x,i). 

Proof. Observe that for any x ^ X and any i E the particles at the sites i_ and i+ 
(call them left and right ones) are consequent. After each iteration of the map T the distance 
between these particles changes by the difference between their velocities, which might take 
values between and v. Since the left particle can be slowed down only by the right one, we 
see that for any moment of time t the distance between the particles can be enlarged at most 
by w\. Indeed, the longest time spent on the same site by the left particle while the right one 
already started moving cannot exceed w. On the other hand, the maximal distance the right 
particle can cover during this time is wv. Thus 

< {i+ + L{x, i+,t)) — (i_ + L{x, t)) < — + wv 

or 

i- — i+ < L{x, i+,t) — L{x, i-,t) < wv. 
Dividing by t and using the definition of the time average velocity we get 

1 1 

\V{x,i,t) -V{x,i+,t)\ <max{-wv, - i_)} *^ 0. (2.9) 

Thus V{x, i) = V{x, Z-|_) for any i. Using the same argument for [i — 1)_ and instead of i one 
extends this result to neighboring particles, and repeating this argument to all particles in the 
configuration. □ 



This shows that if the limit V{x^ i) exists for some site i then the time average velocity is well 
defined, but we still have to check the existence of at least one 'good' site. On the other hand, as 
we have shown the average velocity do not depend on the left tail of a configuration. Therefore 
it is reasonable to introduce for a configuration x G X the notion of its class of equivalence 
X which includes all configurations from X which coincide with x starting from some site, i.e. 
X, y G X if 3fc G : xi = yi^i > k. 



Clearly in the absence of obstacles all particles are moving freely. Therefore to understand 
the dynamics we need to study the motion of 'jams' as the only possible obstacles to the free 
motion of particles. We say that a segment x[n, m] with n < m corresponds to a jam if 

0<Xi<l Vi G {n,... ,m}, (2.10) 

Xn-l < 0, \Xm+l \ > 1. (2.11) 

In other words the jam x[n, m] is a locally maximal collection of consecutive particles having zero 
velocities with the possible exception of the leading one (located at the site m) whose velocity 
is strictly less than 1. For example, in the configuration 

...(a)...(0)...l...(00)l...(000a).l... 

the jams are marked by parentheses. 

The number of particles and their positions in a jam may change with time: leading particles 
are becoming free and some new particles are joining the jam coming from behind. However, 
only one such change at a time might happen, and, in particular, a jam cannot split into several 
new jams. Therefore we can analyze how a given jam changes with time and the main quantity 
of interest for us here is the minimal number of iterations after which the jam will cease to exist. 
Denote by J{t) the segment corresponding to the given jam at the moment t (in this notation 
J(0) is the original jam). Then by the life-time of the jam J we shall mean 

r(J) := sup{t : | J(t)| > 0, t > 0}, (2.12) 

where \A\ is the length of the segment A. 

'Attracting' the preceding particles, a jam plays a role similar to an attractor in dynamical 
systems theory. Therefore it is reasonable to study it in a similar way and to introduce the 
notion of the basin of attraction (notation BA(J)) of the jam J := x[n,m], by which we mean 
the segment x[k,m] with k < n < m of the configuration x containing all particles which will 
eventually at positive time join the jam or can be stopped by a particle from the jam J during 
the time t(J) + l7 Examples of basins of attraction and their dynamics are shown in Fig 01 

Let us introduce the weight W{x[n,m]) of a segment x[n,m] with n < m as 

M^(x[n, rra]) := wNp(x[n,m — 1]) + w{xm) 

= w ■ {m — n) ■ p{x[n, m — 1]) + w{xm)- (2-13) 

The first term in this expression gives the contribution from all particles in the segment except 
for the leading one, which is described by the second term. 

Theorem 2.1 Let v = 1, a < 1 and let J(0) := m] with n < m be a jam. Then its basin 
of attraction BA(x[n, m]) is the minimal segment x[k,rri\, k < n < m for which relations 

Nh{x[k,m]) + l = W{x[k,m]), (2.14) 

Xk-i<l — ia, i = l,2 (2-15) 

^The last assumption excludes the only possibility for a particle to create a new jam when it does not join an 
existing jam J but is stopped at the time r( J) + 1 by the last particle in it which just get the velocity 1 but did 
not start moving yet. Consider an example: .l.a ..01. Here a particle preceding the jam J containing the only 
particle with velocity a is stopped at the next moment of time (and creates a new jam) when the jam J already 
ceased to exist. 
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Figure 4: Dynamics of basins of attraction whose boundaries are marked by square brackets. 

hold true,^ where Nh{x[k,m]) is the number of holes in the segment x[k,m]. Moreover, 

t{J)=W{BA{J)), (2.16) 

i.e. the life-time of the jam is equal to the weight of its basin of attraction, and under the 
dynamics the BA of a jam is transformed to the BA of the remaining part of this jam. 

Relations ()2.14|) . 1)2. 15() give a very simple algorithm to find the left boundary A: of a BA: 
we move to the left until the number of holes in the segment x[k,m] will become equal to the 
weight of this segment minus one. After this we check the condition 1)2. 15() and continue the 
procedure if it does not hold or stop otherwise. 

Proof. First, let us rewrite the relation ()2.14|1 in a more suitable way. From the definition of 
the weight function it follows that for k < m we have 

Np{x[k,m]) = \Wix[k,m])/w]. (2.17) 

On the other hand, we have the trivial identity 

m — k + 1 = Np{x[k, m]) + Nh{x[k, m]), 

which together with 1)2. 17() yields a new relation equivalent to (|2.14l) : 

m-k + 2 = W(x[k,m\) + \W{x[k,m\)/w']. (2.18) 

Consider now in detail the dynamics of a BA of a jam. Two examples are shown in Fig. \^ 
where a = 1/2 and we denote as usual holes by dots and mark particles by their velocities. In 
both examples the jams consist of particles having initially zero velocities, and the right panel 
demonstrates the reason of the inequalities ()2.15|) . which exclude too slow particles that cannot 
join the jam. Observe that if one changes the velocity of the first particle in the right panel 
from to, say a, then the left boundary of the BA goes to the left since the particle immediately 
preceding the left bracket will be stopped by the last particle in the former jam at the moment 
of time when it gets velocity 1. 

Let us show that a basin of attraction of a jam either consists of a single site, or its first site 
is occupied by a hole (i.e. Xk = —1). To do this, consider a function 

(pie) ■.= W{x[m-e,m])+ \W{x[m-e,m])/w'], £€R\. 

According to the definition of the weight of a segment we can write the left boundary of the 
BA(x[n + l,m]) = x[k,m] as k = m — i, where £ is the smallest nonnegative solution to the 
equation 

e + 2 = ^{e). (2.19) 

*When there is no segment satisfying H2.14^ and l|2.15|l we set k = -co, i.e. the BA coincides with the entire 
'left tail'. Note that the relations II2.14II . (12.1511 give only an implicit information about the BA, which differs 
significantly from the result in where the much simpler case a = v has been considered. 



The function is piecewise constant with jumps of amphtude {w + 1) at integer points 
corresponding to the lattice sites i where Xj > (i.e. where there is a particle in the configuration 
x), and (/)(0) = w{xm) + 1- Note that at a jump point the value of cj) corresponds to the upper 
end of the jump. 

There might be two possibilities: 

• Xm G [1 — a, 1). Then = w{xm) + 1 = {w — 1) + 1 = w. Hence the smallest 
nonnegative solution to the equation (|2.19|) occurs at the origin. This case corresponds to 
the BA consisting of a single site. 

• Xm £ [0,1 — a). In this case (/)(0) = w{xm) + l>{w — l) + l = w. Thus the starting point 
^(0) of the graph (j){t) is higher than the starting point of the straight line i + w. Hence 
the first intersection between these graphs occurs at a horizontal piece of the graph of (j){t} 
and thus it corresponds to a hole. 

It remains to take into account the inequalities (|2.15j) . If at an intersection point they are 
satisfied, then the result is proven, otherwise, due to the presence of particles in the segment of 
length 2 immediately preceding to the interval [m — i,m], the graph (f)(£) makes some positive 
jumps of amplitude w + 1 each. Therefore considering the graph after these jumps we are in the 
same position as in the second case considered above and can apply the same argument about 
the intersection point. Checking again the inequalities (|2.15j) and continuing the procedure (if 
necessary) we are finishing this part of the proof. 

Thus if A; < m then the first site k of x[k,m] is occupied by a hole (i.e. Xk < 0), while the 
last site m is always occupied by a particle with the velocity Xm < 1 (since x[n, m] is a jam). 
Therefore after the application of the map T the velocity of the particle at the site m increases 
by a and either this particle becomes free and leaves the jam (if Xm + a, > 1, which yields that 
the new leading particle of the jam is located at the site m' := m — 1), or the velocity remains 
below the threshold 1 and m' := m. In both cases the weight of the segment x[k,m'] decreases 
by 1. Hence, the solution to the equation 1)2. 19(1 with m = m' occurs now at the point i' := £+1 
which yields k' := k + 1, i.e. the left boundary of the BA is shifted by 1 position to the right, 
and the resulting segment x[k' ,m'] is again the BA of the remaining part x[n,m'] of the jam. 

Note also that once we have found the solution to the system (|2.14j) . (|2.15j) then for all 
subsequent time steps the inequalities 1)2. 15(1 will be satisfied automatically, since the particles 
immediately preceding the BA cannot outstrip the left boundary of the BA which is moving 
with the velocity 1. 

Now, since after each iteration of the map T the weight of the BA decreases by 1, we deduce 
that the original weight of the BA is equal to the life-time of the jam. □ 

Let us give also a short explanation to the "physical meaning" of the formula 1)2. 18() . As 
we just have shown after each iteration the left boundary of the BA increases by 1, while the 
right boundary decreases by 1 once per w iterations, and thus the "space-time shape" of the BA 
represents a skewed to the right triangle (see Fig. IIJ. When eventually the BA size vanishes, 
the length of the original BA calculated using these two observations gives the relation ()2.18() . 

This result shows that if all jams in the configuration x ^ X have finite BAs (or finite 
life-times, which is equivalent) then eventually with time all particles will become free. 

3 Validity of the fundamental diagram 

Let us consider the dynamics of jams in more detail. 



Lemma 3.1 Let v = l, a < 1, y £ X and let p+{y) < 71 := (1 + wv) ^. Then there is a 
configuration x G y in which only finite life-time jams may be present. 

Proof. Choose a representative x G y satisfying the assumption that p{x[—n,0]) p+{y) 
and assume contrary to the statement of the Lemma that there exists a jam J = x[n,m] with 
n < m < 00 having an infinite BA. (We do not exclude the case n = — 00 here.) Then by ()2.18|) 
using the same construction as in the proof of Theorem 12.11 for any k < m we get: 

m — k < W{x[k, m]) + [H^(x[A;, m])/w~\ . 

On the other hand, by (O 

Ty(x[/i;, m]) < p{x[k, m]) ■ (m — k + 1) ■ w. 

Hence, 

p{x[k, m]) ■ [m — k + \) ■ w + \p{x[k, rri\) ■ (m — k + > m — k. 
Passing to the limit as k ^ 00, we come to a contradiction: 

P+{x) > {w + l)-^ = ji > p+{x). 

□ 

As we shall see the critical value 71 := {wv + 1)"^ gives the upper bound for the existence 
of a "stable" free-flowing phase. One might think that under the conditions of Lemma 13.11 all 
configurations with densities below the critical value 71 cannot have infinite life-time jams. To 
show that this is not the case consider the configuration y(^) having particles only at sites with 
negative coordinates. Then p{y^^^) = 0, but the segment 00,— 2] of this configuration 

corresponds to a jam, and since it has an infinite number of particles, the life-time of this jam 
is infinite as well. 

Lemma 3.2 Let v = l, a < 1, xEX and let in the configuration x there exist a jam J{t) 
with the infinite BA (i.e., \J{t)\ > Vi > 0). Then each hole in x, located originally to the right 
of J(0), starting from the moment of time when this hole meets with the particle that was the 
leading one in the jam J(0) will start moving periodically by one position to the left exactly 
once per w time steps. 

An immediate extension of this result is the following observation. Let us call the largest 
moment of time at which a given hole meets a particle the life-time of this hole. Then under 
the assumption of Lemma 13.21 the motion of each hole having an infinite life-time is eventually 
periodic with the period w. 

Proof. Observe that since there is at least one infinite life-time jam, then each particle 
belonging to its BA (i.e. located to the left of J) will eventually join the jam J. On the other 
hand, according to the dynamics the leading particle in the jam J is becoming free exactly once 
per w time steps. Thus these free particles are passing holes (i.e. exchanging positions with 
them) exactly once per w time steps (needed for the leading particle in a jam to achieve the 
velocity 1 and to start moving ).9 □ 



®Note that even when v > 1 and some holes may be trapped inside a jam, the situation remains the same due 
to a similar argument. Moreover, whence a particle will join some jam (non necessarily with an infinite BA) and 
will go through it, it will start passing holes exactly once per w time steps. 



This result provides us with the key observation to the calculation of the time-average velocity 
in the jammed phase. 

Assuming that a configuration x G X satisfies the assumptions of Lemma 13.21 let us consider 
in detail the dynamics of a hole located initially at a site j to the right of an infinite life-time 
jam. According to Lemma 13.21 after some finite transient period this hole will start moving 
periodically by one position to the left once per w time steps. Denote by 'dj the duration of the 
transient period and by Kj the number of particles which the hole will meet during this time. 
Then we can characterize the deviation of the movement of the hole from the periodic one by 
the defect of the transient period Dj := Kj — 'dj/w. 

We shall say that a configuration a; G X is ultimately jammed if there is a representative x 
from the same equivalence class in which for any n G there is a jam Jn{t) in the configuration 
X with the infinite BA and Jn(0) > n, i.e. at the moment t = the jam Jn(0) starts further to the 
right from the site n. Since the dynamics of a particle do not depend on particles located to the 
left from it we shall assume (to simplify notation) that in an ultimately jammed configuration 
all sites with nonpositive numbers are occupied by particles with zero velocities. 

Now we are ready to formulate the regularity property mentioned in Introduction. Let in an 
ultimately jammed configuration x holes with nonnegative positions are located at sites {jk}- 
We shall say that the configuration x is regular if the functional 

Reg(x) := limsup \Dj^\/jk (3.1) 

describing statistics of absolute values of normalized defects vanishes. In other words the regu- 
larity means that defects of transient periods can grow only at a sublinear rate. Observe that 
any spatially periodic ultimately jammed configuration is regular^*^ and that the property to be 
ultimately jammed is preserved under the our equivalence relation. 

The following statement describes the limit velocity statistics for high density configurations. 

Lemma 3.3 Let v=l, a<l, xGX and let x be an ultimately jammed configuration.^^ 
Then 

< f— ^ - l) Reg(i). (3.2) 



Proof. Observe first that we are in a position to apply Lemma 13.21 For any z G the total 
distance L{x, i, t) covered by the particle initially located at the site i_ in the configuration 
X £ X during the time i > is equal to the number of holes encountered by the particle 
during this time. For a given t > denote by it the original position (at moment t = 0) in 
the configuration x of the last hole our particle meets during the time t. Then in the segment 
x[i-,it] there are exactly L{x,i,t) holes and, hence, its length it — i- can be found from the 
equation: 

{it - i^) ■ (1 - p{x[i^,it])) = L{x, i, t), 
while the number of particles in this segment is equal to 

Np{x[i^,it]) = {it - i-) ■ p{x[i^,it]). 

Thus 

NMi-^it]) = L{x,i,t) ■ P^'^f-:''^^^ . (3.3) 

1 - p{x[i^,it\) 



^"clearly this property holds if distances between successive infinite hfe-time jams are uniformly bounded. See 
Lemma 13.71 for more general situations. 

^^According to Lemma |3. II if the density is below 71 then all jams have finite BAs. Therefore the assumption 
about the presence of jams with infinite BAs yields the implicit assumption on the density. 



During the time t the hole initiahy located at the site it will meet with all those Np{x[i-,it]) 
particles. Since we are interested in the long term dynamics we may assume that t is so large 
that there is at least one infinite life time jam between the sites z_ and it, which guarantees that 



t>'&i^. Therefore 



w w 



Substituting the obtained relation to equation 13.31 we get 



w * I - p{x[i-,it])' 

which gives 

V{x,i,t) = -L{x,i,t) 



w \p{x[i-,it]) J \p{x[i-,it]) J t 
Observe now that t>it and hence 

|AJA< lAJAt. 

Therefore passing to the upper/lower limit as time goes to infinity and taking into account that 
it *— — * oo we get the result. □ 



Another more intuitive argument for the proof of this statement based on a space-time 
averaging will be discussed in the Appendix. 

The following result demonstrates the importance of the regularity assumption for the con- 
vergence of the limit velocity. 

Lemma 3.4 Let x be an ultimately jammed configuration with Reg(x) > 0. Then even if 
P-{x) = p+{x) = p{x) < 1 the limit points of the time average velocity may differ from the 
value (l/p(x) — l)/w. 

Proof. If the regularity assumption is not satisfied then there exists a sequence of positions of 
holes < ji < ^2 < • • • such that 

lim lAJ/ife = Reg(x) > 0. 

fc— >oo 

Assume contrary to our claim that V{x) = (l/p(x) — l)/w. Then it/t ^II^ V{x) > and 
therefore the second term in the relation (0) describing fluctuations around the time average 
velocity does not vanish with time. We came to the contradiction. □ 



Corollary 3.5 Let x be an ultimately jammed configuration with the density 71 < p{x) < 72 
and with Reg(x) > and let the lower limit of normalized defects differ from the upper one. 
Then using the same argument as above one can show that the values V±(x) differ as well. 

Lemma 3.6 The assumption in Lemma l3 . 31 that the configuration x is ultimately jammed is a 
necessary one and it is always satisfied if 



p-{x) > 72 := (1 + v) \ 



(3.5) 



Proof. If there is only a finite number of jams then particles originating from them might 
meet holes that were never met by other particles before. Assume now that in the configuration 
X only finite life-time jams can be present. Since a free particle should have at least v holes 
immediately ahead of it and since the length of the BA of a jam of n particles exceeds w{n — l), 
we deduce that in any segment of length i there are at most i/(v + 1) particles, provided that 
£ is large enough. 

It remains to show that under the assumption (|3.5|) there is an infinite number of infinite 
life-time jams. Assume that this is not the case, i.e. there exists n G such that in the 
segment x[n,oo] only finite life-time jams are present. Then we can apply the argument in 
the first part of the proof to this segment to demonstrate that this assumption leads to the 
inequality p_(x[n, oo]) < 72. We came to the contradiction since p-{x[n', 00]) = /?_(x[n, 00]) for 
any n' < n. □ 



A more detailed analysis of the regularity functional improves substantially the result of 
Lemma EiSl Moreover, the sufficient condition for the regularity property can be formulated in 
terms of gaps '■= x[nfc,mfc] between successive infinite life-time jams. Define 



Reg(x) := limsup |Gfc|/mfc. (3-6) 

k—*oo 

Lemma 3.7 Let v=l, a < 1, x ^ X and let x be an ultimately jammed configuration. Then 



(a) < Reg(x) < 1/2, 



(b) Reg(a;) > Reg(x), hence Reg(x) = implies Reg(x) = 0, 



(c) Reg(x) = Reg(x) = if p-{x) = p+{x) > 72. 

Proof. Consider a gap between two successive infinite life-time jams Jk,Jk+i- According 
to its definition the gap may contain only free particles and finite life-time jams together with 
their basins of attraction (otherwise those jams will have an infinite life-time). Therefore for 
each particle in the gap there should be at least one hole, which yields the assertion (a). 

To prove the assertion (b) observe that for a hole located initially at a site j we have 
< i^j- < -dj and thus 

\Dj\ = \Kj - ^j/w\ < (1 - l/w)^j. 

Now the statement follows from the observation that for a hole located in the gap Gk the length 
of the transient period cannot exceed \Gk \ + 111. 

Assume now contrary to the last assertion that there is a configuration x (z X with p- (x) > 72 
and Reg(x) > 0. Then there is a sequence of gaps G^ in x such that 



rrik 



Reg(x) > 0. 



According to the proof of the assertion (a) the density of particles in any gap Gk cannot exceed 
1/2. Therefore considering the configuration 2; as a sequence of successive alternative blocks of 
particles Jk and gaps Gk and taking into account that lengths of jams Jk should go to infinity 
with k (otherwise P-{x) < 72 = 1/2 since \Gk\ 00) we see that either Reg(x) = or 

P-{x) < p+{x) (since x consists of alternative linearly growing blocks of two types Jj with 
p{Ji) = 1 and Gi with p{Gi) < 1/2). □ 



From these technical statements we derive the main result about the fundamental diagram. 



Theorem 3.1 Let v = 1, a < 1, x G X. If 

(a) p+{x) < 7i := (1 + wv)^^ then V{x) = v, 



(b) P-{x) > 7i and x is not ultimately jammed (hence P-{x) < 72) then V{x) 

(c) p~{x) > 71 and x is ultimately jammed then 



V±(x) 



1 



1 li; 



-1 



< 



1 



p=p(x) 



1 Reg(x), 



(d) p^{x) = p+{x) = p{x) > 72 then V{x) = - l)w ^ 

Proof. By Lemma l3. II under the condition p+(a;) < 71 only jams with finite life-times may be 
present in x. Consider a partition of by nonintersecting finite BAs of jams in the configuration 
X and their complement. Choose one of those BAs and denote by i the site containing the first 
particle preceding this BA. Even if originally this particle does not have velocity one it will get 
this velocity (and will start moving) at most after w iterations. On the other hand, by the 
definition of the BA this particle (which does not belong to any BA to the right of it) will never 
join a jam after the first moment of time when V{x,i,t) = 1. Finally Lemma 12.21 finishes the 
proof of item (a). 

To prove the existence of the upper branch of the fundamental diagram it is enough to show 
that for any < 7 < 72 := (l + v)~^ there are configurations with particle density 7 and such 
that eventually with time all their particles will become free. To demonstrate this, consider a 
space periodic configuration x £ X with the space period of length 2i in which only even sites 
0, 2, . . . , [2^7] are occupied by particles with velocity 1, while all other sites are occupied by 
holes. Clearly, all particles in this configuration are free and will remain free under dynamics. 

On the other hand, the density of this configuration p{x) = ^ L^^tJ 7- 

The remaining part of Theorem 13 . II describing the lower branch of the fundamental diagram 
follows immediately from Lemmas 13. 3| 13. 6| 13.71 □ 



Obtained results give a complete characterization of the time statistics of a configuration x 
for which the density p{x) is well defined or if it is not well defined but either /0+(x) < 71 or 
p_(x) > 72. Let us show that there are other situations with much wilder behavior. 

Lemma 3.8 Let v = 1, a < 1. Then there exists a configuration y £ X such that P-{y) = 
0, p+{y) = 1 and F_(y) = 0, V+{y) = 1. 

Proof. We shall construct the configuration y as follows. First we set yi := —1 \/i < 0, i.e. we 
fill in the non positive sites by holes. Then we shall fill the remaining sites by alternative blocks 
Bi consisting either only of particles with zero velocities or only of holes. The lengths £j := \Bi\ 
of those blocks we define inductively: 

£1-1, i^:=w'\ h:=w''^'\ ... , ek+i:=w^'=^'% ... 

Denoting := Yli=i S^^ ^k+i = u^"* and Uk+i = + ik+i = + w"'''. Therefore 

p{y[l,n2k+i ) > = ■ — — — > 1. 

n2k+i rik + w'^k 



Using a similar argument, but counting holes instead of particles we get: 



Thus, p{y[l,n2k]) — > 0. This proves the claim about the lower and upper densities. 

Observe now that \/k > the block -Bfc+i corresponds to a jam of length ik+i and that all 
particles to the left of this jam lie in the segment x[l,nfc], whose length is equal to log^^,(^fc+i). 
Therefore by Theorem 12.11 all these particles belong to the basin of attraction of this jam and 
hence will join it with time. On the other hand, since there are no particles at negative sites, 
this shows also that life-times of all jams in the configuration y are finite. Note, however, that 
the life-time of the A;-th jam is of order exp(A;). 

Choose any particle of this configuration and consider how its velocity changes in time. First 
from t = to t = ti >0 the particle might stay in a jam and have zero velocity. Then it is 
becoming the leading one and during w time steps preserves its position but accelerates until it 
get velocity 1. After that the particle starts moving freely with the velocity 1 until it will catch 
up with the next jam. Then it will again stay in a jam having the zero velocity, etc. Due to the 
calculations above the duration of the alternative periods of staying in a jam and free moving 
(interrupted by short periods of acceleration of length w) are growing exponentially. From this 
we get immediately the statement about the lower and upper time average velocities. □ 



4 Stability/instability of the fundamental diagram 

Qualitatively the main difference between the configurations belonging to the upper and lower 
branches of the fundamental diagram is related to the property to be ultimately jammed or not. 
According to Theorem 13. II we need to study the stability of this property only for configurations 
with densities in the region (71, 72) where the two branches of the diagram coexist. The following 
result demonstrate instability of the upper branch of the fundamental diagram, while the stability 
of the lower branch is discussed in Theorem 14.21 

Theorem 4.1 Let v = l, a < 1, xGX and let the density p{x) be well defined and p{x) > 71 
while V{x) = 1 (i.e. the configuration x belongs to the upper branch of the fundamental 
diagram). Then there exists an ultimately jammed configuration y (z X which differs from x on 
a set S having zero density. 

Proof. We shall construct the perturbed configuration y as follows. Set £ := and for for all 
k G for which x^ > consider a sequence of numbers S^/ ■= p{x[(. + l,fc]). We denote by 
ni the value of k which gives the first local maximum of Sk/ for which Sk/ > 71 • Then we set 
H. := ni and continue the procedure to find the number 722 giving the local maximum exceeding 
the value 71, etc. (see Fig. El- Eventually we shall have a monotonically growing sequence 
ni 00, from which we can further choose a subsequence fii satisfying the assumption that 
\hi — nj+il > 2*. Observe that the set of integers {ni}i has zero density. 

Let yj := for all j <= and for positive j we set yj = Xj except for the sites hi where 
we set yfn ■= 0, i.e. we have changed the velocities of the particles at the sites hi to for all i. 
Since n, is the point of a local maximum, the site hi is the leading point of some jam in y. Our 
claim is that the BA of this jam covers the entire segment [hi-i + l,hi]. Assume that this is not 
the case, i.e. there exists N < hi such that the BA corresponds to the segment y[hi — N + l,hi] 
and denote by M the number of particles in this segment. Due to the local monotonicity of the 
sequence Sk/ in the segment preceding the point hi the strict inequality 71 < M/N holds true. 
On the other hand, by Theorem 12. II we have 

N = W{y[hi -N + 1, hi]) + [W{y[hi -N + l,hi])/w\ = wM + M. 

Thus M/N = = 71, which contradicts to our assumption. 



71 
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Figure 5: Averages along a configuration 



Now, since the density of differing sites is equal to zero (due to the choice of hi) while the 
configuration y satisfies the assumptions of Lemma we get the result. □ 



Clearly, changing back the zero density set of sites in the ultimately jammed configuration y 
constructed in the proof above we get the configuration x in which all particles will eventually 
move freely. Therefore one might expect that all ultimately jammed configurations with densities 
in the interval (71, 72) are also unstable. The following statement shows that this is not the case. 

Theorem 4.2 Let v = 1, a < 1. There exists a regular configuration x ^ X with 71 < p{x) < 72 
and constants < yl < i? < 00 such that if a configuration y G X differs from x at most at 
A sites in any segment of length B then y is also regular. In other words, there is an open 
neighborhood of the configuration x consisting only of configurations corresponding to the lower 
branch of the fundamental diagram. 

Proof. Choose positive integers N, M large enough such that 71 < M/ [N + M) < 72 and 
consider a spatially periodic configuration x whose main period x[l,N + M] consists of N 
consecutive holes and M consecutive particles having zero velocities. Then the density p = j^jj 
of this configuration is well defined and this configuration evidently satisfies the conditions of 
Lemma 13.31 

Our aim is to show that for a small enough (but still nondegenerate) perturbation the 
property to have a countable number of infinite life-time jams arbitrary far to the right is 
preserved. Moreover we shall see that there exists a finite A'^i such that any segment of length 
A^i in the perturbed configuration y intersects with some infinite life-time jam. Clearly additional 
particles in y (compared to x) should not worry us and it is enough to consider only the case when 
under the perturbation some particles are removed from the configuration x or their velocities 
are changed. 

Choose n = k{N + M) with k E Z-|- large enough and assume on the contrary that all jams 
in y[l,n] have finite life-times and their BAs lie completely in this segment (i.e. they do not 
include the site 0). Denote by mo the number of free particles in the segment ?/[!, n] and by mi 
the number of particles belonging to jams in this segment. 

By Theorem 12.11 we know that if / is a finite segment corresponding to a basin of attraction 
then W{I) > wNp{I) and Nh{I) >W{I). Thus 

|/| = Npil) + Nh{I) >iw + l)NpiI). 



Taking into account that each free particle is followed by a hole we conclude that 



2mo + {w + l)mi < n. 



On the other hand, 
Thus 



niQ + mi = Np{y[l, n]) < np. 
mo ^ iw + l)p - 1 



> 



> 0. 



n w — 1 

Observe now that in order to get a free particle in the configuration y we need to change at least 
two sites in the configuration x. From this one immediately can find the constants A, B setting 



A ■-- 



{w + l)p - 1 
w — 1 



■ B 



and choosing B := 2k{N + M) to be so large that A>1. 



□ 



Using the same argument one can prove stability for a regular configuration for which instead 
of space periodicity one assumes that the gaps between infinite life-time jams are uniformly 
bounded. 

Let us show now that a single site perturbation of a configuration consisting only of free 
particles can create a jam whose length grows linearly with time. For v=l, a = 1/2 consider 
a configuration x having free particles (with velocity 1) at all even sites and holes at all others. 
We perturb this configuration only at the origin setting the velocity at site to zero and thus 
creating a jam of length 1. According to Theorem 12.11 the BA of this jam is infinite. On the 
other hand, after each iteration a new particle joins this jam from the left while only once per 
two iterations the leading particle leaves the jam, which proves its linear growth. 

A partial result in the direction of a measure-theoretic interpretation of the notion of regu- 
larity for configurations in the density region (71,72) gives the following Lemma. 

Lemma 4.1 Let a configuration x has a density < p{x) < 1. Decompose x into alternative 
blocks of particles Ai and holes Bi: x = . . . A1B1A2B2 ■ ■ ■■ Then 

\An \ + \-Bn\ n— >oo „ 

^=i{\A\ + m) ~" • 

In other words, the lengths of blocks can grow only at a sublinear speed. 

Proof. Let us prove that lengths of blocks of holes cannot grow at a linear speed. Assume 
contrary to our claim that there exists a constant < C < 1 and a sequence of integers < 
ni < mi < n2 < m2 < . . . with < ^mi for all i such that all sites in the segments x[ni,mi] 
are filled by holes while the segments x[mi, rij+i] are filled by particles. Since the density of the 
configuration x is well defined we have: 

p{x[l,k]) p{x). 

On the other hand, 

^(^[1,^^]) = ^^P(^i^^^^]) < Cp{x[l,nS (P{X), 



which contradicts to the assumption that the density of x is well defined. A similar argument 
applies to blocks of particles as well. □ 



Recalling the connection between the functionals Reg and Reg we see that the set of regular 
configurations do not differ much from the set of all configurations having densities even in the 
"hysteresis" region and we expect that for any nontrivial translationally invariant measure on 
X the set of regular configurations is a set of full measure. 



5 Appendix: Proof of Lemma 13.31 via space-time averaging 

Lemma 13.31 can be proven also in a more intuitive way using a kind of a space-time averaging. 
Since this approach explains the connection between the time and space averaging we shall 
discuss it additionally to the formal proof given in Section |31 

Observe that by the definition of the map T in the case v = 1 we explicitly have the mass 
conservation law, i.e. after each move one particle exchanges its position with one hole. Therefore 
we can give a kind of 'physical derivation' of the relation ()3.2() . Assume that the (time) average 
velocity V{x) is well defined and choose a segment of the configuration x of length £, provided 
^ ^ 1. Then the number of particles in this segment is equal to Np := l{p{x) + o{\ / €)) , while the 
number of holes is equal to := 1{1 — p{x) + o{l / 1}) . We are in a position to apply Lemma 
hence the holes are moving with the average velocity to the left. Thus in total Np particles 
will move to the distance tNpV{x) during the time t, while in total A'^^ holes will move to the 
distance tN^/w during this time. Due to the mass conservation these quantities coincide and 
thus 

U{p{x) + o{l/l))V{x) = U{1 - p{x) + o{l/l))/w. 
Dividing by t and l and passing to the limit as ^ ^ oo, we get 

1 - p{x) 



V{x) 



wp{x) 



To make this argument rigorous, we need some additional work to be done. We start with 
the space-periodic case and restrict ourselves to just one period. For any given time interval the 
total shift of all particles in the 'spatial period' is well defined. Dividing this total shift by the 
time and by the number of particles in the 'spatial period' and passing to the limit as time goes 
to infinity (which exists due to the behavior of holes) we get the 'space-time' average velocity. 
On the other hand, since V{x, i) does not depend on i we deduce that our 'space-time' average 
is, in fact, just V{x). Let us prove the last statement. Denote the period length by For any 
moment of time t we have the following relation: 

^ t 
V l>o{xi) ■ L{x, i, t) = £{1 - p{x)) ■ - + R{t), 
^-^ ~ w 

i=l 

Here (due to the usage of the indicator function) the summation is taken over all particles in the 
'spatial period', the remaining term R{t) cannot exceed l/w on absolute value, and the value 
£{1 — p{x)) is equal to the number of holes in the spatial period. Dividing by we get 

' - 1 

>o(xO • V[x, i, t) = t{l - p{x)) ■ - + 0(l/t), (5.1) 



El 



■ 1 ^ 



where 0{\/t) means a term of order 1/t as t — > oo. 

On the other hand, by the inequality 1)2. 9j) for any two particles initially located at sites i,j 
of the 'spatial period' and any moment of time t we have: 



1 1,. .2 



\V{x,i,t)-V{x,j,t)\ <£max{-u;v, <r/t. 



I 1 ^ I I I > t 

hit') hit') 

Figure 6: Dynamics of traffic jams at moments of time t and t' > t. 



provided that £ is large enough. Thus 



(x) • Vix,i.,t) - ^l>o(xi) ■ Vix,i,t) 



i=l 



< i^t, 



where the value ipix) is equal to the number of particles in the 'spatial period'. 

Therefore passing to the limit as i ^ cxd in the relation ()5.1() and using the inequality above, 
we get 

t^oc ipyx) w \pix) J 

which finishes the proof in the space-periodic case. 

In the general non space-periodic case one also can follow a similar argument. According to 
our assumption in the configuration x there are finite jams with infinite life-times, denoted 
by rectangles in Fig. |B1 Here the parameter t indicates the moment of time. The segment which 
starts immediately after the site occupied by the leading particle of the jam Jiit) and finishes 
at the site occupied by the leading particle of the jam Jj+i, we denote by Then one can 

proceed as follows. 

(a) The number of holes in hit) is invariant with respect to the dynamics. Indeed, if a hole 
was located at time t to the right of a jam J(t), then for any moment of time t' > t when Jit') 
still exists this hole should be still to the right of Jit'). The reason is that the right boundary 
of a jam changes only when its leading particle becomes free. However, the new leading particle 
of the jam still remains to the left of our hole. 

(b) The number of particles in can change due to dynamics at most by 1 in both sides. 
Indeed, the leading particle of each the jams Jj(t) is becoming free periodically with the period 
w, which is the only way how the number of particles in liit) can change. 

According to these properties one can think about the holes inside of each segment as 
a kind of a pump which is pushing the particles through itself with the constant 'productivity' 
equal to 

Nhilijt)) l_ 
Npiliit)) ' w ' 

which leads to the desired result. 



Concluding remarks 

A deterministic generalization of the Nagel-Schreckenberg traffic fiow model (as well as of the 
slow-to-start model) with the real valued acceleration has been proposed and studied analytically. 
It has been shown that macroscopic physical properties of the model under study such as the 
time average velocity of particles in the flow depend crucially on the density of particles in 
the flow and these results are described in terms of the corresponding fundamental diagrams. 
Moreover, we have shown that instead of space (or space-time) average velocities used earlier 
in traffic models the time average velocity can be considered. It is proven that this quantity 



characterizing the behavior of an individual particle coincides for all particles in the flow and 
in distinction to the space average velocity (which wildly fluctuates in time) leads to a proper 
statistical description of the dynamics. 

One of the problems for future analysis is the question what happens when the initial config- 
uration is random or the dynamics is weakly perturbed in random or deterministic sense (e.g., 
random traffic lights are taken into account). The main idea in the paper is that without pertur- 
bations one can predict the asymptotic behavior of an individual configuration (and even of an 
individual particle in it). Clearly integrating the results according to some initial distribution 
one gets predictions for a family of random initial configurations. One would expect that such 
questions might be of interest only in the 'hysteresis' region where both branches of the funda- 
mental diagram are present. Even there according to Theorem 4.1 the influence of the upper 
branch should be negligible. Another point is that it might be possible that under arbitrary 
small perturbations the upper branch of the fundamental diagram will disappear^^ and that the 
non-regularity of initial configurations should not matter in the random setting. However to 
check these predictions one needs to study a probabilistic version of the model. 

Finally let us mention that it is important to distinguish between peculiarities related to a 
trivial divergence of various series due to the fact that we deal with the infinite phase space 
(infinite lattice) and true finite-size phenomena taking place even in the case of spatially 
periodic configurations (or a finite lattice). The latter issue seems more important and our 
result proves that the hysteresis phenomenon exists even in this case. 

The author is grateful to anonymous referees for very helpful comments. 



^^From |1U| it follows that if v = 1, o = 1/2 and random perturbations are applied only to jams but not to free 
particles the upper branch exists at least for spatially periodic configurations. 
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